Effects of the vortices and impurities on the nuclear spin relaxation rate in iron-based 

superconductors 
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The effects of magnetic vortices and nonmagnetic impurities on the low energy quasiparticle 
excitations and the spin-lattice relaxation rate are examined in the iron-based superconductors for 
the s±-, s- and d-wave pairing symmetries, respectively. The main effect of the vortices is to enhance 
the quasiparticle excitations and the spin-lattice relaxation rate for all symmetries, and leads to a 
T'^ dependence of the relaxation rate followed by a nearly T-linearity at lower temperatures. This 
enhancement can only be seen for the s±- and d-wave symmetries in the presence of nonmagnetic 
impurities. These results suggest that the s±-wave and d-wave pairing states behave similarly 
in response to the magnetic field and nonmagnetic impurities, therefore it may be impossible to 
distinguish them on the basis of the measurements of spin-lattice relaxation rates when a magnetic 
field and/or impurity scatterings are present. 
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I. INTRODUCTION 



Recently much attention has been payed to the newly 
discovered iron arsenide superconductors^^- which dis- 
play superconducting transition temperature as high as 
more than 50K, appear to share a number of gen- 
eral features with high- Tc cuprates, including the lay- 
ered structure and proximity to a magnetically ordered 
statejii^ii These observations suggest that the conven- 
tional phonon-mediated pairing mechanism appears to 
be unlikely and the magnetic correlations may be rele- 
vant for superconductivity. So far, unconventional super- 
conductivity with pairing symmetry s± mediated by the 
interband spin fluctuations has been proposed by a num- 
ber of theories for this layered iron superconductors.—"— 
Although, such a popular proposal can explain some ex- 
perimental findings, the situation is complicated by the 
power-law temperature dependence of the spin-lattice re- 
laxation rate T^^ ^ T" below Tc with the doping de- 
pendent n varying from 6 to l.SiJ^"^^ Because a simple 
theoretical analysis shows that the relaxation rate T^^ 
should exhibit an exponential temperature dependence 
for a fully gaped superconductors, while power law rela- 
tion Tf ^ - holds in the presence of line nodes in the 
SC gap. Moreover, several experiments have reported 
the evidence of a residual density-of-state (DOS) at zero 
energy in the SC state, where the temperature depen- 
dent relaxation rate deviates from the Tj~ ^ ~ relation 
and exhibits T," 



ture 
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' T behavior at very low tempera- 
The disparity between the theoretical pro- 
posal and the experimental fact presents a puzzle that 
should be resolved for the determination of the SC pair- 
ing symmetry in these high-Tc superconductors. 

The relaxation rate has been studied in the pres- 
ence of the impurity-enhanced quasiparticle scattering 
at zero field in the band representation, and it shows 



the s± = Aq cos(fca;) cos(fcy)-wave SC pairing 
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though the enhancement of the quasiparticle scattering 
has been well included in the previous studies, the re- 
lationship between the SC gap and the impurities ad- 
dressed in the band representation were debated contro- 
versiallyi2L2^ More importantly, actual spin-lattice relax- 
ation rate measured in NMR experiments are conducted 
under a magnetic field of several Tesla, where the relax- 
ation rate shows the temperature dependence T^^ 
and even Tf ^ ~ j^j^^^^ ^t the same time, H. J. 
Grafe et al. has found that the superconducting vortices 
contribute to the spin-lattice relaxation rate when the 
magnetic field is perpendicular to the conducting plane 
but not for the parallel directionJ^i^ In fact, it is un- 
clear up to now whether the different value of n in the 
temperature dependence of T^^ is caused by a change of 
the pairing symmetry with doping, a disorder scattering 
effect in an s± gap superconductor, or other unknown 
mechanism. In view of these theoretical and experimen- 
tal facts, it is necessary to compare and contrast different 
pairing symmetries by studying the effects of impurities 
and magnetic field on the spin-lattice relaxation rate in 
the analysis of the standardized procedure extracting the 
gap symmetry. 

The purpose of this study is to present such a con- 
trastive study on three different superconducting (SC) 
pairing symmetries with s±-, s- and d-wave, respectively. 
For this purpose, we calculate the DOS and the spin- 
lattice relaxation rate in the Fe-based superconductors 
in the presence of magnetic field and nonmagnetic im- 
purities by self-consistently solving the Bogoliubov-de 
Gennes equations based on the simple two-orbital model. 
It is shown that the characteristic low energy quasi- 
particle excitations depend on the gap functions in the 
presence of magnetic vortices and nonmagnetic impuri- 
ties. The magnetic vortices contribute significantly to 
the spin-lattice relaxation rate and lead to the relation 



the power law temperature dependence Tj^ ^ ~ T'^ for T-^ ^ ^ y3 followed by a nearly linear dependence at lower 
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temperatures for all symmetries. While in the presence 
of nonmagnetic impurities, this enhancement of Tj"^ can 
only be seen for the s±- and d-wave symmetries. Thus, 
in the presence of magnetic field and nonmagnetic impu- 
rities, the s±-wave pairing behaves much like what the 
d-wave does. 

The remainder of the paper is organized as follows. In 
Sec. II, we introduce the model Hamiltonian and carry 
out analytical calculations. In Sec. Ill, we present nu- 
merical calculations and discuss the results. In Sec. IV, 
we make a conclusion. 



orbitals \ti \ = 1 as the energy unit, we have. 



II. THEORY AND METHOD 

We start with an effective two-orbital model that takes 
only the iron dxz and dyz orbitals into account)^ By as- 
suming an effective attraction that causes the supercon- 
ducting pairing and including the possible interactions 
between the two orbitals' electrons, one can construct an 
effective model to study the vortex and impurity physics 
of the iron-based superconductors in the superconducting 
state: 



H — Hq + Hpair ■ 



The first term is a tight-binding model 



(1) 



Ha 



ij.ap.cr 



which describes the electron effective hoppings between 
sites i and j of the Fe ions on the square lattice, includ- 
ing the intra- {tij^aa) and inter-orbital {tij^a,i3,0L ^ /3) 
hoppings with the subscripts a, j3 (a, = 1,2 for xz 
and yz orbital, respectively) denoting the orbitals and a 
the spin, cj creates an a orbital electron with spin 
a at the site i {i = {ix,iy)), and fi is the chemical po- 
tential. The magnetic field is introduced through the 
Peierls phase factor e"^*^ with ipij — JJ^' ■A-(r) ■ dr in 

the vortex state, where A — {—Hy, 0,0) stands for the 
vector potential in the Landau gauge and $o = hc/2e 
is the superconducting flux quantum. In the case of the 
SC state with impurities, we randomly select the half 
number of total sites, at where the random disorder po- 
tentials Ui uniformly distributed over [~U,U] are set in. 
Then U is a parameter to characterize the strength of 
the disorder. The hopping integrals are chosen as to cap- 
ture the essence of the density function theory (DFT) 
results.'^^ Taking the hopping integral between the dyz 
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(3) 

Here, x and y denote the unit vector along the x and y 
direction, respectively. 

The second term accounts for the superconducting 
pairing. Considering that a main purpose here is to ad- 
dress the effects of the magnetic vortices and nonmag- 
netic impurities on the spin-lattice relaxation rate in the 
iron-based superconductors, we take a phenomenological 
form for the intra-orbital pairing interaction. 



Hpair 



cl„^c].+h.c.) (4) 



with Vij as the strengths of effective attractions. 

Thus, we obtain the Bogoliubov-de Gennes equations 
for this model Hamiltonian^^ 
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a,<T (^"/3,a); ^la.a («"/3,ff) are the Bogohubov quasi- 
particle amplitudes on the j-th site with corresponding 
eigenvalues En- 

The pairing amplitude and electron densities are ob- 
tained through the following self-consistent equations, 

n 

*""^(2S^^) 

"^c.t = E l<a,tl'/(^n) 

71 

n 

The site-averaged DOS N{E) is calculated by 

^(^)= -^EE[i<c.,tiv'(i?n-i?) 

i n.a 

+ K^A^f\En + E)l (8) 
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where, / (E) is the derivative of the Fermi-Dirac dis- 
tribution function with respect to energy. The nuclear 
spin-lattice relaxation rate is given by-** 

R{rt,ri,) = Imx+ -(ri,rj/,il7„ n + ii])/{n/T)\Q^o 

XTrTf'{Er,)S{E„-E„,). (9) 

Here, = u^^ + and V," = + vlp. We choose 
Yi = Til because the nuclear spin-lattice relaxation at a 
local site is dominant. Then the site-dependent relax- 
ation time is given by Ti(r) = 1/R{r,r) and the bulk 
relaxation time Ti = {l/N) i]^ri(r). 

In numerical calculations, Vij is chosen to give a short 
coherence length of a few lattice spacing in the SC state 
being consistent with experiments.'^'^ Under the condi- 
tions Vij ^ 2.0, /i = 1.2 at temperature T = 1 x 10^'"', 
the filling factor n — ^ „{na,a) / {NxNy) — 1.9 and the 
coherent peak of the SC order parameter in the DOS is 
at Amaa; ^ 0.25. Thus, we estimate the coherence length 
^0 ~ EFa/\/S.max\ ~ So^^ with a being the Fe-Fe dis- 
tance on the square lattice. Due to this short coherence 
length, presumably the system will be a type-II supercon- 
ductor. To study the vortex states, we employ the mag- 
netic unit cell with size iV^; x Ny = 48 x 24 that accommo- 
dates two magnetic vortices, unless otherwise specified. 
In view of these parameters, we estimate the upper crit- 
ical field Bc2 ~ lOOT. Therefore, the model calculation 
is particularly suitable for the iron-based type-II super- 
conductors such as CaFei_a;Coa;AsF, Euo.7Nao.3Fe2As2 
and FeTei-ajS^;, where the typical coherence length 
deduced from the experiments is of a few lattice spacing 
and the upper critical field achieves as high as dozens of 
Tesla^ 



III. RESULTS AND DISCUSSION 

Since no final consensus on the SC pairing symme- 
try has yet been achieved, we choose three possible sin- 
glet pairing symmetries, i.e., the most popular sign- 
reversed s±-wave, the on-site s-wave and the d-wave 
symmetries with their respective gap functions As_,_ — 
Aq cos(A;^) cos(fcy), A^ = Aq and A^ = Ao[cos(fca;) — 
cos(fcy)] to carry out the contrasting study. At the end 
of this section, we will also touch on another possibility 
of the s±-wave with angular variation along the electron 
Fermi pockets, which will be referred to as s^-wave. In 
order to obtain comparable values of critical temperature 
Tc in the self-consistent calculations, we set respectively 
Vij = 1.6 for s-|--wave, Vij = 2.0 for s-wave, and Vij = 1.8 
for d-wave pairing symmetries. 

To begin with, we briefly summarize the site-averaged 
DOS spectra N{E) in the uniform SC state at T = 1 x 
10^^ as shown by the solid lines in Fig. 1. For the 
s±- and s-wave symmetries, no node exists in the gap 
along the Fermi surfaces. Correspondingly, the full gap 



structures can be seen in the DOS as shown in Figs. 1(a) 
and 1(b). In the d-wave symmetry, the gap function has 
line nodes at — ky, which cross the hole pocket but do 
not intersect the electron pocket. As a result, the DOS 
consists of a small V-shaped gap structure at very low 
energy and a U-shaped gap structure at higher energy as 
shown in Fig. 1(c), which exhibits a difference from that 
in high- Tc cuprates. 

Next, we show the DOS in the magnetic vortex state 
and in the presence of nonmagnetic impurities as pre- 
sented by the dotted and dashed lines in Fig. 1, respec- 
tively. In the vortex state, the application of a magnetic 
field will induce the quasiparticle flow around the vortex 
core, such that the nodal line will appear in the otherwise 
fully gaped SC state due to an additional Doppler shift in 
the quasiparticle energy, giving rise to the V-shaped gap 
structures for the s±- and s-wave symmetries, as shown 
in Figs. 1(a) and 1(b). In contrast, in the presence of im- 
purities, the DOS shows the V-shaped structure for the 
s-|--wave symmetry but the U-shaped for the s-wave sym- 
metry [Figs. 1(a) and 1(b)]. This is due to the fact that 
the impurity potential has intra- and inter-band compo- 
nents for the multiband materials. The intra-band com- 
ponents scatter fermions that have the same sign for the 
s-|--wave SC order parameter and therefore do not affect 
the superconductivity. Whereas the inter-band compo- 
nents scatter fermions with opposite SC order parame- 
ters, thus have the pairing breaking effect. As a result, 
they yield an obvious decrease in and simultaneously 
introduce the V-shaped feature in DOS. In the case of 
the s-wave symmetry, no sign change in the SC order pa- 
rameters occurs on both the electron and hole pockets, 
so that no obvious pairing breaking effect is induced by 
the impurities and the U-shaped DOS is untouched. 

For the d-wave symmetry, the small V-shaped super- 
conducting gap is filled by the low energy quasiparticles 
induced by either vortices or impurities, resulting in a 
pseudogap-like U-shaped feature with finite density at 
Fermi energy, as indicated by the dotted and dashed lines 
in Fig. 1(c). The line nodes existing on the hole-pocket 
in the d-wave symmetry make it vulnerable to the im- 
purities, resulting in the disappearance of the small V- 
shaped structure. However, the full gap opening on the 
electron-pocket is robust against impurities, so that the 
pseudogap-like U-shaped feature is obtained. In the vor- 
tex state, since the quasiparticles induced by the vortex 
come preferably from the line nodes, the pseudogap-like 
U-shaped structure remains. The results in Fig. 1 indi- 
cate that the characteristic low energy quasiparticle ex- 
citations depend on the gap functions in the presence of 
magnetic vortices and nonmagnetic impurities. 

Now, we turn to the discussion of the temperature (T) 
dependence of the nuclear relaxation rate. For the uni- 
form SC state, both the s±- and s-wave symmetries pro- 
duce a power law relation T^^ ~ below Tc and it 
evolves into an exponential dependence at very low tem- 
perature, as shown by the solid lines in Figs. 2 and 3)22 
which are the consequence of the full-gap DOS in Fig. 1. 
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In the vortex state, due to the similar V-shaped DOS for 
both the s±- and s-wave symmetries, T^^ changes it's T 
dependence to below Tc and becomes nearly propor- 
tional to T for the s-t-wave symmetry while T^'^ for the 
s-wave symmetry at low temperature, as denoted by the 
dotted lines in Figs. 2(a) and 3(a), respectively. [The 
dotted lines in Figs. 2(a)-5(a) show the results for mag- 
netic unit cell with size 48 x 24, while the dash-dotted 
lines in these figures the results for magnetic unit cell 
with size 40 x 20.] 

We note that both the T'^ dependence and the low-T 
feature with nearly linear slop are reminiscent of the ex- 
perimental observations r^ii^ where the T-linear depen- 
dence has been regarded as the evidence for a residual 
density of states at zero energy in SC state. As men- 
tioned above, although the s±-wave pairing is basically 
the full gap, the application of a magnetic field will cause 
the quasiparticle flow around the vortex core, such that 
the nodal line will appear in the SC state due to an ad- 
ditional Doppler shift quasiparticle energy, giving rise to 
the V-shaped DOS and dependent relaxation rate. In 
view of this, we may expect that the slop of tempera- 
ture dependent T^^ should be insensitive to the strength 
of the magnetic field, which has already been observed 
in experiment 1^ This is evident by the comparison be- 
tween the dotted- and dash-dotted-curves shown in Fig. 
2(a), where the magnetic field for the dash-dotted-curve 
is about 1.5 times as large as that for the dotted-curve. 

A striking difference between the s-t-wave and the s- 
wave states in the T-dependence oiT^^ appears when the 
effect of the impurity scattering is considered, as shown 
in Figs. 2(b) and 3(b). For the s±-wave state, Tj~^ devi- 
ates gradually from to a overall behavior at weak 
disorder such as C/ = 1 , then to a T linear dependence at 
low temperatures as the disorder strength U is increased 
to about U — 1.5, which depicts the sensitivity of the 
SC order to the impurities and is in accordance with the 
DOS results. On the other hand, Tf^ changes little upon 
the introduction of the disorder for the s-wave symmetry. 
We also notice that Tc is reduced substantially for the 
s-|--wave symmetry at the moderate and even the weak 
disorder strength, as can be seen from the lower-T shift of 
the inflexion point on the curves in Fig. 2(b). This again 
reflects the fact that the sign-reversed s±-wave pairing 
is fragile against the nonmagnetic impurities, which has 
also been predicted theoretically by adopting a more so- 
phisticated orbital model and conformed experimentally 
in the spcciflc heat and resistivity measurements 

In Fig. 4, we present the results for the d-wave sym- 
metry. Unlike the formers, the overall T-dependence of 
T^^ in the uniform SC state roughly follows rela- 
tion due to the existing of nodal line in the gap struc- 
ture. When a magnetic field is applied, T^^ exhibits a 
T-linear dependence at low temperatures, though a 
behavior following Tc still remains. We notice that this 
trend is rather robust against the magnitude of the mag- 
netic field, as the results for two cases are nearly the same 
though their magnitude differs about 1.5 times, as shown 



in Fig. 4(a). Thus, for the three different pairing states, 
T^^ exhibits nearly similar T-dependence in the mag- 
netic vortex state, due to the presence of quasiparticles 
induced by the vortex. In the presence of impurities, T^^ 
for both the d-wave and the s±-wave symmetries exhibit 
a consecutive change from the to T-linear dependence 
when temperature is decreased, which contrasts with that 
for the s-wave state. 

We notice that the actual multiple Fermi surface sheets 
in iron-based superconductors are not exactly reproduced 
by the two-orbital model, so the general structure of the 
gap in the s±-wave channel may involve the angular vari- 
ation along the electron Fermi pockets and have the form 
A/,. = Ao{cos(fcj:) cos{ky) + h[cos{kx) + cos(fcy)]}, where 
the factor h measures the strength of the angular depen- 
dent variations along the electron Fermi pockets.''*'"'*'' In 
the case oi h > 1, there will be accidental nodes along 
the electron Fermi pockets,^^ which will be focused here. 
Such accidental nodes are refiected in the V-shaped DOS 
in Fig. 1(d) and nearly T"* dependence of Tj~^ for the 
typical results with /i = 1.5 in the uniform SC state [Fig. 
5]. This result differs from the dependence for the 
symmetry imposed nodal behavior such as in the case of 
the d-wave state. Although the T-dependence of T^^ in 
the vortex state is much like that of the s±-wave case, it 
is less influenced by the impurities [see Fig. 5]. 

IV. CONCLUSION 

In conclusion, we have investigated the effect of the 
magnetic field and nonmagnetic impurities on the DOS 
and the spin-lattice relaxation rate T^^ in the iron-based 
superconductors. It is shown that the characteristic site- 
averaged DOS depends on the gap functions in the pres- 
ence of magnetic vortices and nonmagnetic impurities. 
The magnetic vortices have a significant contribution to 
the spin-lattice relaxation rate and lead to the relation 
Tf 1 - T3 followed by a nearly T-linear dependence at 
low temperatures for all three symmetries (s±-, s- and d- 
wave) considered here, though in the clean uniform state 
a dependence for the d-wave symmetry differentiates 
from the others with a T^ dependence. In the presence 
of nonmagnetic impurities, this enhancement oiT^^ can 
only be seen for the s±- and d-wave symmetries, whereas 
it is almost unaffected for the s-wave symmetry. Our 
results suggest that it is impossible to distinguish the 
s±- and d-wave symmetries on the basis of the measure- 
ments of spin-lattice relaxation rates when a magnetic 
field and/or impurity scatterings are present. 
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FIG. 2: T-dcpcndcncc of T^^ shown in the double logarithmic 
chart for s±-wavc symmetry in the vortex state (a), and in 
the presence of nonmagnetic impurities (b). The dotted and 
dash-dotted curves denote different strength of the magnetic 
field (see text) and the impurity scattering. The results for 
the uniform SC state are also plotted with solid line in each 
figure. 



FIG. 1: Site-averaged DOS for the s±-wave in (a), the s- 
wave in (b), the d-wave in (c), and the s^-wave in (d) (see 
text). The DOS in the uniform SC state is plotted with solid 
lines. The results in the magnetic vortex state and those in 
the presence of nonmagnetic impurities are shown with the 
dotted and dashed lines, respectively. 
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